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ABSTRACT 

We present a study of the dust properties of 12 classical and superbubble HII regions in the 
Large Magellanic Cloud. We use infrared photometry from Spitzer (8, 24, 70, and 160 /zm bands), 
obtained as part of the Surveying the Agents of a Galaxy's Evolution (SAGE) program, along with 
archival spectroscopic classifications of the ionizing stars to examine the role of stellar sources 
on dust heating and processing. Our infrared observations show surprisingly little correlation 
between the emission properties of the dust and the effective temperatures or bolometric magni- 
tudes of stars in the HII regions, suggesting that the HII region evolutionary timescale is not on 
the order of the dust processing timescale. We find that the infrared emission of superbubbles 
and classical HII regions shows little diff'erentiation between the two classes, despite the signifi- 
cant differences in age and morphology. We do detect a correlation of the 24 ^m emission from 
hot dust with the ratio of 70 to 160 fim flux. This correlation can be modeled as a trend in the 
temperature of a minority hot dust component, while a majority of the dust remains significantly 
cooler. 

Subject headings: HII Regions — dust - ISM: bubbles 



1. Introduction 

The life-cycle of interstellar dust remains a 
largely under-constrained process. Until recently, 
attempts to characterize the evolution of dust 
have been impeded by uncertainties regarding 
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dust's composition, physical state, emission pro- 
cesses, and illuminating radiation field. All of 
these parameters must be known prior to draw- 
ing any observationally-based conclusion on dust 
evolution, and before recent results from the 
Spitzer Space Telescope emerged, these parame- 
ters had been largely uncertain as well. A consen- 
sus is gradually forming around some of these 
properties, with the modeling of the "uniden- 
tified infrared bands" as vibrational emission 
from polycycli c aromatic hydrocarbons (PAHs; 



Leger fc Pugetl Il984t lAUamandola et al.l llQSSi) . 



and models of the size distribution and cross sec- 
tions of the PAHs and sma ll grains are now better 
defined (|Li fc Draindl200l[ ). 

Many Spitzer programs have taken advantage of 
these new dust models to provide insight on many 
different classes of objects. Previous works have 
char acterized the dust content of entir e galaxies 
fe.g. lDale et al.ll2005l: lOraine et al.ll2007l) and indi- 



vidual star forming regi ons (e.g. ICompiegne et al 
I2OO8 : Berne et al. 20071 among many others), but 
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consistent analyses of larger samples of HII regions 
have be en somewhat lacking. One exception is the 
work of Gordon et al. ( 2008h . which used a sam- 
ple of HII regions in MlOl to probe metallicity- 
dependent dust processing. Their findings linked 
the ionization parameter of HII regions with a 
decrease in the PAH emission feature equivalent 
widths, which is evidence for the destruction of 
PAHs. 

Our work seeks to relate the properties of the 
stars in HII regions to their dust content and emis- 
sion. The radiation output from these massive, 
luminous stars dominates the heating of the dust, 
and the stellar SED and luminosity may affect 
the dust emission. Evolved massive stars are also 
thought to be an important source for dust pro- 
duction, and stellar wind shocks and supernovae 
can also be important for processing dust. Our 
work leverages a previously compiled sample of 
HII regions with known spectral types for the ion- 
izing stars. The sample includes regions with ages 
spanning from 1-5 Myr, including both classical 
HII regions and evolved, multi-SN superbubbles. 
This range of objects enables us to estimate the 
dust processing timescales. We seek to determine 
which of these varied parameters, if any, dominate 
the dust properties and emission. In Section [2l 
we describe the observations and data reduction; 
Section |3] contains our interpretation and model- 
ing of the observed infrared ratios; Section |4] dis- 
cusses the relation between stellar properties and 
infrared emission; and we explain the applicability 
of our work to other observations in Section [S] 

2. SAGE Observations 

Our sa mple of HII regions was drawn from the 
objects in lOev &: Kennicutt ( 1997 ). The sample is 



listed in Table[l] with identifiers from lDavies et al 



( 1976L ; hereafter DEM). The sample contains 8 
superbubbles and 4 classical HII regions, and all of 
the objects in this sample have spectroscopically- 
determined spectral types for their ionizing stars. 
In Table [1] we have listed the spectral type of the 
hottest star in the association, which dominates 
the stellar luminosity and spectral energy distri- 
bution in the region. We have also included the 
total stellar bolometric magnitude for the ionizing 
stars, determined by summing the luminosities of 
all stars spectroscopically observed by the work 



referenced in the last column. The bolometric cor- 
rections had been applied by the authors of each 
cited work except for DEM 1 99, for which we used 
the bolometric corrections of Martins et al. 
for O and B stars and ICrowtheil (|2007t ) for Wolf- 
Rayet stars. With the inclusion of superbubbles, 
our sample spans a wider range of object ages, 
which could potentially demonstrate dust evolu- 
tion. Many of the superbubbles also have super- 
novae and shocks which could drive dust process- 
ing. 

The infrared data were obtained from the 
Spitzer Legacy program. Sur veying the Agents of a 
Galaxy's Evolution (SAGE: [Meixner et al.ll2006[ ). 
The program observed an area of ~ 7° x 7°, cov- 
ering the LMC in aU IRAC (3.6, 4.5, 5.8, and 
8^m) and MIPS (24, 70, and 160 fim) bandpasses. 
The individual observations were calibrated, com- 
bined, and mosaicked by the SAGE team, with a 
pipeline similar to that o f the GLIMPSE sur vey 
(see Meixner et al. 2006: Benjamin et al.|[2003l for 
further details). 

2.1. Photometry 

In order to compare surface brightnesses of the 
target objects at different wavelengths it is nec- 
essary to convolve the images to the lowest res- 
olution of the observations, which in our case is 
the 160 fim band (full- width at half-maximum of 
40"). Due to diffractive features in the instru- 
ment's point-spread function (PSF), such as Airy 
rings and spikes, a convolution with a Gaussian 
kernel would not produce sufficiently accurate re- 
sults. Instead we used th e con volution kernels 
described in iGordon et akl (|2008l ). which are cal- 
culated from model PSFs. This convolution also 
alleviates the need for aperture corrections when 
computing band ratios, since the aperture correc- 
tion will be a constant factor for all bands. 

We followed the prescription of iHelou et al 



(j2004l ) for using the 3.6 /im images to remove the 
stellar flux from the 8 and 24 /im images. We com- 
puted the non-stellar flux (F"*^) from the observed 
flux (-Fi/) in each band using 

F;'"(7.9^m) = F^(7.9/im)-0.232F^(3.6/im) (1) 

and 

F^"(24/im) = F^{24: fim)~0m2F^{3.6 ^lm) , (2) 
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where the coefficients are from lHelou et al.l (|2004l ). 
These constants were determined by comparing 
the mid-infrared emission of stellar population 
synthesis models to their 3.6 /itm emission. The 
infrared emission of these models is not strongly 
sensitive to star formation history or metallicity 
and should be generally applicable, but they are 
still approximate. Using the 3.6 /im images for 
this subtraction also assumes that the 3.6 fim band 
includes no dust emission, even though a small 
amount of dust emission is often present in this 
band. A visual check of the subtraction confirmed 
the efficacy of the subtraction at removing numer- 
ous faint point-sources from the images. Since the 
contribution from stellar sources at 71 and 160 /im 
is negligible, we performed no correction on those 
bands. When discussing the 8 /.tm and 24 ^m pho- 
tometry in what follows, we will be referring to 
these corrected non-stellar fluxes. 

We performed aperture photometry on the ob- 
jects, with apertures defined by hand. In all 
cases, we set the apertures such that they in- 
cluded all of the Ha flux that could reason- 
ably be associated with the target object. The 
Magellanic Cloud Emission-li ne Survey (MCELS; 
ISmith fc MCELS Teamlll998l ) Ha imaging of the 
LMC was used for this comparison. Figure[T]shows 
infrared and Ha images of all of the superbubbles 
in our sample, and Figure [2] shows the classical HII 
regions. In both figures the photometric apertures 
are outlined in white. 

2.2. Background and Error Estimation 

Our photometry is potentially contaminated by 
both diffuse emission from the Milky Way and by 
emission from unrelated dust in the LMC in the 
foreground or background of our objects (hereafter 
collectively referred to as "background" emission) . 
We subtracted this background emission from our 
objects, measuring the background level with an 
annulus spanning 14 to 17 arcminutes from the 
center of each region. We used the mode (the most 
commonly occurring value after binning the values 
into 0.01 MJy/steradian bins) of the pixels in this 
annulus as the background level. The mode was 
chosen rather than the mean in order to reduce the 
sensitivity of the background estimation to bright 
structures from neighboring clouds. We note that 
it is impossible to completely account for emission 
unrelated to our objects that may appear within 



our target apertures, since much of the dust emis- 
sion varies on size scales smaller than our aper- 
tures. Though this should not cause a systematic 
error in our photometry, it will still necessitate a 
level of caution in interpreting values for specific 
target objects. 

Given that it is possible for this background 
contamination to exist despite our efforts to re- 
move it, we would like to know the extent to which 
it may compromise our measurements. We there- 
fore need to know the variance we would mea- 
sure with our apertures if they were placed ran- 
domly and not on our target objects. By placing 
a number of "error estimation" apertures of the 
same size as our object apertures and measuring 
the variance in integrated light amongst them, we 
can thus obtain some estimate of the uncertainty 
in our photometry. This still does not allow us 
to measure the background's contribution to any 
single object's photometric measurement any bet- 
ter than we did with the background annulus, but 
it does inform us of the uncertainty in that mea- 
surement. The flux densities for all of our objects, 
corrected for stellar and background emission, are 
listed in Table These flux densities have not 
had an aperture correction applied, so they are 
likely to be systematically low. This does not 
affect the band ratios used in this work though, 
since the necessary correction factors are identical 
for all bands as a result of convolving each image 
to a common PSF. Contamination by unremoved 
background is the dominant source of uncertainty 
in our photometry and the error bars on the fig- 
ures and in Table [2] reflect that. 

3. Dust Emission and Modeling 

To allow for comparisons between objects with 
different luminosities, we constructed three differ- 
ent flux rati os from the d ata, following the meth- 
ods of .Draine fc Lil (|2007t ). These are 



7.9 — 



i^F;^^{7.9tim) 



and 



vF^(n ^m) + vF„{im^i\n) ' 
vF^{n fini) + i^F^{160fiin)' 
vF^{71 ^m) 



R7 



i^F^{160 iim)' 



(3) 
(4) 

(5) 
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P7.9 is primarily a tracer of the aromatic emis- 
sion commonly associated with PAHs, while P24 
tends to trace hot, thermalized dust. Both of 
these are normalized by the total thermal dust 
emission, traced by the quantity i/F^(71 fim) + 
vF^{160 ^m). The ratio Rn measures the temper- 
ature of the dust that has reached thermal equilib- 
rium, which is generally comprised of larger grains 
with radii greater than 100 A. Since the temper- 
ature at which these grains equilibrate is deter- 
mined solely by the amount of radiation they re- 
ceive, the 70/160 ^m ratio is also useful as a tracer 
of the total radiation content of the region. 

Figure [3] shows the effect of radiation inten- 
sity on both the PAH and VSG populations. The 
top panel shows that Py.g is largely uncorrelated 
with the total radiation content of the region. 
This is expected as the PAH emission is the re- 
sult of stochastic heating of grains: the level of 
the PAH emission is directly proportional to the 
radiation intensity, but the overall PAH emis- 
sion spectral sh ape does not vary with the radi- 
ation intensity / Praine fc Li Since the ra- 
tio P7.9 is normalized by the far-infrared emission 
{vF^(ll /im) -f vF,j(\QQ /im)) which is also propor- 
tional to the radiation content of the region, the 
ratio is therefore largely independent of the inten- 
sity of the radiation. 

The lower panel in Figure [3] indicates that P24 
ratio increases with radiation intensity. Since the 
8 /im data show that PAH emission is largely un- 
correlated with the 70/160 /im ratio, it is clear that 
the trend in 24 /im must be the result of thermal- 
ized dust, rather than the PAH emission within 
the 24 /im bandpass. At low temperatures, the 
thermal dust emission is detected only in the 70 
and 160 /xm bands, but for dust exposed to high ra- 
diation intensities, the thermal peak can reach into 
the 24 /im band. We therefore hypothesize that as 
the overall luminosity of the region increases, the 
peak radiation intensity seen by the dust and con- 
sequently the temperature of the dust also increase 
in a correlated way. 

Which parameter determines the position of an 
HII region along this 24 /im sequence, and which 
parameters determine the locus of the sequence 
itself? Since we hypothesize that this trend is 
driven by the effects of the illuminating radiation 
field, and not a change in dust composition, we ex- 
amine parameterizations of the radiation that can 



best reproduce this correlation. 

Our simple models of the sequence in 24 /im 
emission are shown in Figure S) In this figure 
we have plotted P24 -- OA4P7.9, which was sug- 
gested bv lDraine et al.l ( 2007 ) as a way to remove 
the contribution of PAHs to the 24 /im band, and 
leave only the flux from thermalized grains. This 
minimizes any effects caused by variations in PAH 
abundances. Figure 2] shows the same tight corre- 
lation in this measure of the dust emission as that 
for P24 alone, uncorrected for PAH emission. 

Using absorptio n cross-sectio n s and grain size 
distributions from Draine &: Lil (2007), we con- 
structed a variety of models to compare to our 
observed sequence. Our models compute only the 
emission from the thermalized dust. In thermal 
equilibrium the dust grain emission must equal ab- 
sorption, hence 

[°° 47r /■°° 

/ Cabs (a, A)uA d\= — / Cabs (a, \)Bx (T) dX, 
Jo c Jo 

(6) 

where Cabs (a, A) is the absorption cross-section as 
a function of grain radius a and wavelength A, u\ 
is the radiation e nergy density, and B\{T) is the 
Planck func ti on (iLi fc Drainel lioOll ) . We follow 
Draine fc Li ( 20071 ) in parameterizing the radia- 
tion field as ux — Uu^^^^^, where we use U 
as a scaling coefficient on the mean interstella r 
radiation field u^^^^^ from iMathis et al.l (|l983l ). 
This parameterization was chosen primarily to aid 
comparison to other works; since the grain is as- 
sumed to reach an equilibrium temperature, only 
the energy content of the radiation field affects the 
dust emission and not the spectral energy distribu- 
tion. For a given radiation field, we compute the 
temperature T for each grain size such that the 
equation of equilibrium (Equation [S]) is satisfied, 
then use that temperature to compute the infrared 
emission of the dust grain. We then sum this 
emission with the emission from gr ains of other 
sizes, using the size distribution from lLi fc Draiii3 
()200lh . 

In the simplest model, all of the dust in the HII 
region is heated by a single radiation intensity. We 
then computed a set of models spanning a range of 
radiation intensities {U), shown by the red, dotted 
line in Figure SI This simple model shows that as 
the parameter U increases, the model's position in 
the P71 — P24 plot moves towards the upper right. 
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This trend follows a similar slope to that of our ob- 
served set of HII regions, but the models exhibit a 
70/160 /im ratio that is much higher than those of 
the observed regions. This suggests that though 
some dust must be exposed to high radiation in- 
tensities to reproduce the 24 /im emission, cold 
dust is also needed to reproduce the observed low 
70/160 /im ratios. These dust components could 
be physically separate and have distinct tempera- 
tures, or they could form a continuum spanning a 
range of temperatures, but dust must be present 
at both high and low temperatures. 

We therefore extend our models to include a 
mass fraction 7 of dust exposed to high radi- 
ation intensities, and a mass fraction (1 — 7) 
that is exposed to a fixed low radiation inten- 
sity [/min- We have a d opted this parameteriza- 
tion from iDraine fc Lil (j2007n . but note that we 
do not include the power law distribution of radi- 
ation intensities that was used in their work; each 
of the dust components sees only a single radia- 
tion intensity. Namely, the dust is heated by two 
radiation intensities, with (1 — 7) C/min of the to- 
tal emission arising from the cold dust (heated by 
Umin) and jU of the total emission arising from 
the hot dust (heated by U). We ran a grid of 
models over this parameter space (7, J/mm) for a 
range of U and sought to minimize the chi-square 
between the models and the observations. The 
best-fitting model had the parameters 7 = 0.0015 
and J/niin — 0.4, shown by the blue, dashed line in 
Figure [4] and with tick marks indicating the posi- 
tion oi U ~ 200 and U — 800 along the sequence. 
The small value of 7 required to fit the observa- 
tions suggests that the bulk of the dust inside our 
apertures is relatively cold, potentially originating 
from the photodissociation region surrounding the 
ionized gas. This cold dust is at temperatures of 
roughly 15 K, but grains of different sizes equili- 
brate at different temperatures and thus there will 
be a range of dust temperatures. For comparison, 
grains exposed to U — 600 equilibrate between 
roughly 40 and 55 K. 

It is important to remember that these models 
are merely suggestive, and cannot be considered a 
definitive explanation of the observed trend in P24 • 
We have here sought to explain the observations 
in terms of a single parameter model, where the 
only difference between objects is the maximum 
radiation intensity seen by a small fraction of the 



dust. It is worth noting that although we do not 
consider the derived parameters (7, Umim U) to 
precisely represent the actual radiation intensities, 
the U values inferred from this simple model (i.e. 
U > 100) are indeed consistent with the premise 
that the 24 /xm emission is mainly from dust at- 
taining equilibrium tempe ratures. A s shown in 



Figure 13 and Figure 15 of iDraine fc Lil (|2007l ). at 
U > 100 the 24 fim emission is dominated by dust 
with equilibrium temperatures. The fact that such 
a simple model is able to reproduce the observed 
trend is certainly interesting, but should not be 
taken as a conclusive determination of the radia- 
tion field illuminating these objects. Our model is 
merely useful for illustrating that the observed se- 
quence is set by the hottest dust component, and 
the temperature that this dust reaches determines 
the HII region's position along the sequence. Sim- 
ilarly, it is beyond the power of our data to suggest 
exactly what amount of dust is in the hot or cold 
components of the HII region itself since we could 
not definitively remove unrelated background or 
foreground emission from our measurements. 

4. EfTects of Stellar Properties 

Figure [5] shows our three infrared ratios plotted 
against the spectral type of the hottest ionizing 
star and the total bolometric magnitude of all the 
stars in the region. The top two panels show again 
that there is little systematic change in the PAH 
emission between objects. The center and lower 
panels can be used to diagnose the dust heating. 
What the data show is that neither the peak tem- 
perature nor the average temperature of the dust 
changes significantly as a function of stellar spec- 
tral type or total region luminosity. 

This is not to say that the total dust emission 
is completely unaffected by the stellar luminosity; 
since they are both ratios they do not probe the 
total dust luminosity. Instead the 70/160 /zm ratio 
tells us that in regions with higher stellar luminosi- 
ties, the dust does not have systematically higher 
temperatures. Looking at the effect of the bolo- 
metric luminosity on the 24 /im emission in the 
middle panel of Figure \5\ we see that the stellar 
luminosity also does not influence the peak tem- 
perature of the dust. In Table[T]we have computed 
^m/^boi, which is the ratio of the measured in- 
frared emission (assumming 24, 70, and 160 /xm 
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bandpasses, but without correcting for emission 
outside these bands) to the bolometric luminosity 
of the stars in the region. Most of the regions have 
values of LiR/iboi between 5 and 15%, though this 
is a lower bound and integrating over the entire in- 
frared spectrum will yield higher values. The fairly 
constant value supports our picture that the total 
infrared emission necessarily increases to accom- 
modate increased stellar emission but, as shown 
above, the dust temperature does not appear to 
be strongly affected. 

We can estimate the temperature difference 
that would be expected from the range of bolo- 
metric luminosities spanned by our objects, using 
basic radiative equilibrium (Equation [5]), which 
implies that the stellar luminosity and dust tem- 
perature are related as L^, oc T^^^^, if the dust 
is like a black body. But because the far in- 
frared dust absorption cross-section is approxi- 
mately Cabs K A~^, the dust temperature is actu- 
ally related to the stellar luminosity by L^, oc T^^^^ ■ 
Our objects span roughly a factor of 100 in lumi- 
nosity, which corresponds to a factor of 2 increase 
in dust temperature. We therefore would only ex- 
pect a modest trend in 70/160 /im, which would 
be further obscured by the large intrinsic scatter 
between objects. Note that this calculation also 
assumes that the dust in all objects is located the 
same distance from the illuminating stars. If the 
dust in more luminous regions was systematically 
closer or further from the stars, then the relative 
dilution of the radiation field may also effect emis- 
sion. 

5. Discussion 

It is important to highlight the fact that the 
correlation between the 70/160 /im ratio and 
24 /xm emission we have observed exists despite 
the numerous parameters of the HII regions that 
could affect their dust content and emission. Ob- 
jects in our sample cover a wide range of mor- 
phologies and illuminating stellar sources. This 
could generate variations in the dust temperature 
distribution and resulting infrared emission, but 
this is not observed. Instead, the objects sit on a 
well-defined sequence with only minimal modifi- 
cation by geometric effects. 

By including in our sample both classical HII 
regions and more evolved superbubbles, our sam- 



ple also spans a range of ages. The superbub- 
bles are generally older HII regions that have had 
time to evacuate a shell of material, and it would 
be reasonable to expect the dust in these ob- 
jects to have undergone more processing than their 
younger counterparts. Additionally, superbubbles 
are often old enough for supernovae to have oc- 
curred within them, which could affect the dust 
content. We do not observe any evidence of such 
effects, though, since in both 8 /xm and 24 /xm the 
two classes of objects exhibit no clear differentia- 
tion from each other. This could be a consequence 
of our significant photometric uncertainties, but is 
compatible with grain processing being completed 
on timescales shorter than the typical age of our 
classical HII regions (1-2 Myr), or on timescales 
longer than the typical age of the superbubbles (4 
- 5 Myr). 

Our results can also explain observations 
of HII regions in other galaxies. Galaxies in 
the Spitzer Infrared Nearby Galaxies Survey 
(SINGS) have exhibited a decreased 8/24 /xm ra- 
tio in clumps corresponding to large HII regions, 
while the 8/160 Atm ratio remains fairly smooth 
(|Bendo et al.ll2008l) . The 8/24 /im ratio is roughly 
a factor of two less in the clumps than it is in the 
rest of the galaxy, a factor which can be easily ac- 
counted for by our 24 /im observations. Our mod- 
els therefore imply that the observed clumping in 
8/24 /im can be explained by higher dust temper- 
atures in these regions, without having to invoke 
variation in dust abund ance. Ano t her s tudy of 
SINGS galaxies by .Cal zetti et al] (|2003) found 



that in the integrated light of galaxies there is 
correlation between 24 /tm and Paschen a emis- 
sion, which suggests that the 24 /im emission is 
indeed related to star formation. Similarly, star- 
burst galaxies also show enhanced 24 fj,m emission 
that is attributed to hotter dust (e.g., Hanish et 
al. 2010). 

Since our work deals only with the integrated 
light from the target objects, we are not able to 
distinguish between emission from the ionized re- 
gion and emission from the photodissociation re- 
gion (PDR). Emission from PAHs has typically 
been fo und within the PD R and not the ionized 
region (Kas sis et al. 2006*) , and consequently the 



properties of the HII region and the PAH emission 
may not be closely linked. We are also unable to 
say if the cold dust emission we observed is from 
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the ionized region or the PDR. With its reduced 
radiation intensity, the PDR would be the logi- 
cal source for the cold dust emission, but emission 
from the ionized region could be present as well. 



Gordon et alj (|2008l ) find a decrease in the PAH 



equivalent widths for HII regions as a function of 
their "ionization index" , which is dominated by 
the ratio of [Ne III]/[Ne II] emission lines. The ion- 
ization index is directly linked to the temperature 
of the ionized gas, which could be responsible for 
the dust processing. We should see a similar effect 
in comparing our superbubbles, which have low 
ionization parameters due to geometric dilution of 
the radiation field, with our classical HII regions. 
But as described above, we find no systematic dif- 
ferences in 8 /xm emission between the two classes 
of objects. Our objects span a com parable range o f 
ionization parameters to the iGordon et al. ( 2008h 
objects, and it is thus surprising that all of our 
objects exhibit so little variation in P7.9. 

6. Conclusions 

We find that the correlation between the 
70/160 /im ratio and 24 /im emission is seen in 
both HII regions and superbubbles. We have qual- 
itatively reproduced this trend with a simplified 
dust model that interprets this effect purely as 
the result of variations in the dust temperature be- 
tween objects. That is, increasing an object's dust 
temperature alters its position in the P24 — i?7i di- 
agram (Figure H]) along the trend we observe. We 
therefore cannot infer grain processing or variation 
in grain composition. 

Additionally, we do not detect any correlation 
between the dust emission ratios and the spectral 
type of the ionizing stars in the regions or the to- 
tal luminosity of all stars in the region. There 
is also no differentiation in infrared emission be- 
tween classical and superbubble regions, despite 
the significant differences in age, morphology, and 
ionization. This suggests that we are not observ- 
ing processing or composition differences on the 
evolutionary timescale for superbubbles. 

This research is based on observations made 
with the Spitzer Space Telescope, which is oper- 
ated by the Jet Propulsion Laboratory, Califor- 
nia Institute of Technology under contract with 
NASA. Support for this work was provided by 
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NSF grant AST-0806476. We would like to thank 
K. Sellgren and A. Witt for helpful discussions, 
along with the anonymous referee for their insight- 
ful comments. 
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Table 1: Stellar Parameters for Target Objects 





Superbubble 


Hottest ionizing 


Total bolometric 




Lib./ Lboi 


Reference 


Name 


or classical 


star 


magnitude 








DEM 25 


S 


09 V 


-8.88 


10 


0.058 


Oev (1996) 


DEM 31 


S 


WN6 


-11.6 


24 


0.011 


Oev d996) 


DEM 34 


C 


03 III 


-13.1 


49 


0.029 


Parker et al. ('1992) 


DEM 50 


S 


06.5 V 


-10.7 


12 


0.073 


Oev (1996) 


DEM 106 


S 


06.5 V 


-11.4 


16 


0.055 


Oev d996) 


DEM 152 


S 


05 III 


-12.4 


64 


0.087 


Oev & Massev ('1995) 


DEM 192 


S 


WC5 


-12.3 


53 


0.0071 


Oev & Smedlev ('1998) 


DEM 199 


C 


WR 


-11.8 


35 


0.15 


Garmanv et al. C1994) 


DEM 226 


s 


06.5 V 


-10.5 


10 


0.091 


Oev (1996) 


DEM 243 


c 


09 V 


-11.0 


17 


0.091 


Oev (1996) 


DEM 301 


s 


03 1 


-10.8 


24 


0.055 


Oev (1996) 


DEM 323 


c 


03 III 


-12.2 


27 


0.13 


Massev et al. fl989) 



"The number of stars whose luminosities are summed to compute the total bolometric magnitude. 



Table 2: Measured Infrared Fluxes for Target Objects 







IRAC4 


MIPS24 


MIPS70 


MIPS160 






8.0 Aim 


24 ^m 


70 /im 


160 fim 


Name 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


DEM 


25 


2.2 ± 4.5 


4.3 ± 2.5 


87 ± 69 


220 ± 200 


DEM 


31 


2.9 ± 1.9 


3.4 ± 1.4 


230 ± 39 


520 ± 89 


DEM 


34 


1.4 ± 16 


1.6 ± 13 


35 ± 360 


90 ± 1000 


DEM 


50 


28 ± 4.8 


28 ± 2.3 


690 ± 48 


1300 ± 130 


DEM 


106 


45 ± 2.3 


110 ± 2.0 


870 ± 42 


1600 ± 100 


DEM 


152 


150 ± 8.6 


530 ± 9.9 


3900 ± 160 


4300 ± 350 


DEM 


192 


6.0 ± 1.4 


23 ± 1.2 


320 ± 25 


390 ± 63 


DEM 


199 


130 ± 24 


440 ± 52 


3800 ± 660 


5300 ± 1100 


DEM 


226 


24 ± 15 


37 ± 18 


760 ± 300 


1200 ± 600 


DEM 


243 


37 ± 6.8 


84 ± 3.6 


1100 ± 79 


1900 ± 230 


DEM 


301 


20 ± 3.3 


30 ± 1.8 


480 ± 57 


1100 ± 140 


DEM 


323 


120 ± 100 


230 ± 100 


3800 ± 1600 


9700 ± 4300 
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Fig. 1. — Images of the superbubbles in our sample in 8, 24, 70 /im and Ha, with our object apertures 
overlaid. The scale bar in t he lower left of each linage is one arcminute, corresponding to 15 pc at the 
distance of the LMC (52 kpc, ISzewczvk et~alll2008l ). 
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Fig. 2. — Same as Figure [TJ but showing the classical HII regions in the sample. 
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8.0/im 
(70/im + 160/im) 




0.1 1 10 

70/160/im 



Fig. 3. — Top: PAH emission as a function of radiation intensity of the region, showing httle correlation. 
Blue squares are superbubbles, while red circles are classical HII regions. Bottom: 24 //m emission shows a 
clear increase as the radiation content of the region increases. Symbols are as above, and the solid line is a 
fit to all points. 
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Fig. 4. — The 24 thermal emission component shown against the 70/160 ^m ratio. Superbubbles are 
blue squares, and classical HII regions are red circles. The best-fit dust model is shown by the blue dashed 
line, which contains a small fraction of hot dust determined by 7. The red dotted line shows a model where 
all of the dust is hot 7 = 1, and the green solid line shows a model with only a miniscule fraction of hot dust 
(7 = 10~^). Tick marks are shown on the best-fit model at [/ = 200 and U — 800 to give an idea of the 
range of radiation intensities (C/) needed to span the observations. 
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Fig. 5. — The infrared emission ratios are plotted against spectral type of the hottest ionizing star in each 
region (left), and against the total bolometric magnitude of massive stars in the region (right). Superbubbles 
are marked with blue squares, while non-superbubbles are marked with red circles. Objects that contain 
Wolf-Rayet stars are circled. Slight horizontal offsets have been added to points in the left panels for clarity. 



